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We derive Bardeen-like expressions for the transmission probabilities between two multi-probe 
mesoscopic conductors coupled by a weak tunneling contact. We emphasize especially the dual role 
of a weak coupling contact as a current source and sink and analyze the magnetic field symmetry. 
In the limit of a point-like tunneling contact the transmission probability becomes a product of 
local, partial density of states of the two mesoscopic conductors. We present expressions for the 
partial density of states in terms of functional derivatives of the scattering matrix with respect to 
the local potential and in terms of wave functions. We discuss voltage measurements and resistance 
measurements in the transport state of conductors. We illustrate the theory for the simple case 
of a scatterer in an otherwise perfect wire. In particular, we investigate the development of the 
Hall-resistance as measured with weak coupling probes. 
PACS numbers: 61.16 Ch, 72.10 Fk, 73.20 At 



I. INTRODUCTION 



Tunneling from a small metallic tip or from a suitable 
mesoscopic contact into a sample is a powerful means for 
the structural analysis of surfaces on an atomic length 
scale Q . In the typical arrangement a two-terminal mea- 
surement is performed using the tip (contact) as a source 
and the sample as a current sink. Modeling the tip as a 
spherical symmetric object (s-wave) and using Bardeen's 
approach Q the current flowing from the tip into the 
sample was found to be proportional to the local density 
of states (LDOS) of the surface at the position of the tip 
[^,|| . In this article we investigate arrangements in which 
the sample is so small that the phase coherence length 
exceeds all sample dimensions. The sample is connected 
to several contacts, which makes it possible to investigate 
it in a transport state. In particular, multiterminal resis- 
tance measurements become possible. In one interesting 
configuration two of the contacts of the sample act as 
current source and sink and the STM (or contact) serves 
as a voltage probe (see Fig. |^). 




FIG. 1. Probing the voltage with an STM tip on a meso- 
scopic wire. Wire and tip are coupled via a single bond con- 
necting point Xq in the wire to point Xq in the tip with a 
coupling element t. 

Pioneering experiments using such a configuration 
have been undertaken by Muralt et al. Q and Kirtley 
et al. P . Theoretically weak coupling contacts to small 



conductors were considered already by Engquist and An- 
derson to find the resistance of a scatterer. The dis- 
cussion of Engquist and Anderson led to Landauer's re- 
sult 1^ which expresses the resistance of a scatterer in 
terms of the ratio of its reflection probability R and its 
the transmission probability T. Landauer's derivation 
makes no appeal to measurement probes, but following 
Engquist and Anderson subsequent multichannel gener- 
alizations of this result by Azbel and Biittiker et al. 
JlOt were also given interpretations as weak coupling mea- 
surements by Imry . It was later pointed out that the 
discussion of Engquist and Anderson is not exact but 
neglects Friedel-like oscillations, generated by the inter- 
ference of incident waves reflected at the scatterer [ p^ . 
Interest in the multichannel generalization of Azbel faded 
when it was noticed that it could not account for the 
magnetic field asymmetry observed in metallic diffusive 
wires |l3j. Furthermore, in mesoscopic physics experi- 
ments probes typically are massive and cannot be treated 
as a weak coupling measurement. A general multitermi- 
nal approach which treats all probes on an equal footing 
was introduced by Biittiker |l4| . This approach, which 
expresses resistances as rational functions of transmission 
probabilities, also permits a discussion of weak coupling 
probes. But to our knowledge a detailed investigation of 
weak coupling probes based on this general approach has 
not been carried out. Below we will present such an anal- 
ysis in which the weak coupling probes and the sample 
are described with one overall scattering matrix. 

While the initial experiments by Muralt et al. and 
Kirtley et al. |^ are difficult to interpret improvements 
in sample preparation techniques and in low temperature 
STM will hopefully lead to a resumption of such measure- 
ments. In two terminal configurations Friedel oscillations 
in the equilibrium electron density near steps | p^ , p^ and 
in atomic chorals [0 have been observed and there is 
no principal reason why a similar resolution could not be 
achieved in the measurement of the transport state. A 
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detailed discussion of potential oscillations near a barrier 
has been given by Chu and Sorbello |18|. These authors 
also find a marked difference between the potential mea- 
sured at the weak coupling probe and the electrostatic 
potential in the bulk of the sample. 

It is the purpose of this work to derive general weak 
coupling formulae for multiterminal conductors. Of par- 
ticular interest is a formulation in which the dual role 
of the weak coupling contact, which can act either as a 
current injector (or sink) or as a voltage probe, appears 
in a manifestly reciprocal way. We start from the multi- 
terminal formulation of Ref. (ijjl^ which is valid for all 
types of contacts. 




FIG. 2. Two mesoscopic multiterminal conductors (called 
system A and system B) are weakly coupled by a coupling ma- 
trix Tab which connects a region of conductor A to a region 
of conductor B. 

The most general arrangement which we consider is 
depicted in Fig. ||. Two conductors A and B are weakly 
coupled by a coupling matrix Tab- Each conductor is 
connected via ideal leads to several electron reservoirs. 
For the special case that the two conductors are only 
coupled by a single tunneling path from to Xq like it 
is the case in Fig. |^ we find that the transmission prob- 
ability from contact a of conductor A through the weak 
coupling contact to contact (3 of conductor B is given by 
a generalized Bardeen-like expression, 



a 



(1) 



Here i is a coupling energy, va{xq ^ a) is the injectivity of 
contact a into point x^ and VBiPjXQ) is the emissivity 
of point x^ into contact /3. The injectivities and emissiv- 
itics po| , ^ arc only a portion of the LDOS and will be 
denoted as local partial density of states (LPDOS). Be- 
low we give analytical expressions for these densities of 
states in terms of functional derivatives of the scattering 
matrix of conductor A and B and in terms of wavefunc- 
tions 1^,^ . Here we mention only that in the presence 
of a magnetic field B the injectivity and emissivity in 
each conductor are related by a reciprocity relation: The 
injectivity from a contact a into a point a; in a magnetic 
field B is equal to the emissivity of this point into contact 
a if the field is reversed. 



v{B; x, a) = v{—B] a, x) 



(2) 



Consequently the transmission probability given by Eq. 
(||) manifestly obeys the Onsager-Casimir symmetry 



Tai3{B) = Tpa{-B). Eq. d) is one of the central re- 
sults of this work. It can be used to find the resistances 
in an arbitrary multiterminal geometry with the help of 
formulae that express these resistances as rational func- 
tions of the transmission probabilities. 

Like transmission probabilities the LPDOS are eval- 
uated in the equilibrium electrostatic potential. They 
are thus affected by interaction (screening, etc.) only 
to the extend that interactions determine the equilib- 
rium electrostatic potential. On the other hand the LP- 
DOS are not sufficient to determine the actual electron 
distribution in the sample. The change in the charge 
density in response to an increase of the electrochem- 
ical potential, say at contact a by d/io,, is given by 
dn{x) = i'{x, a)dfia + dni„dix) where driindix) is the in- 
duced charge density gen erated by the non-equilibrium 
electrostatic potential [po|] . 

Below we introduce the Hamiltonian formulation 
p3|-p5| of the scattering matrix and express the LPDOS 
using this approach. The LPDOS are also expressed in 
terms of scattering states. We present the derivation of 
Eq. (|l|) starting from the full scattering matrix of the 
weakly coupled system. We then discuss a number of 
applications and the relation of our results to the earlier 
work mentioned above. The magnetic field dependence 
of our resistance formula is illustrated on a ballistic con- 
ductor with a barrier. 



II. HAMILTONIAN FORMULATION OF THE 
LOCAL PARTIAL DENSITIES OF STATES 

In a first step we derive expressions for the local par- 
tial density of states in terms of the Hamiltonian of the 
sample. Let us consider a mesoscopic conductor which 
is connected via ideal leads to electron reservoirs. We 
assume that, at the Fermi energy, we have in each lead 
a, Na open channels. The Hamiltonian approach p^-p^] 
starts with a formal division of the Hilbert space into two 
parts, the open leads and the compact sample region. 
The Hamiltonian can then be written as 



n = K + H + W + W'^ 



Here 



N No, 



K ~ \am){am\EF 

a— 1 rn—1 

is the Hamiltonian of the isolated leads. 



H 



J2\x){x'\H, 



(3) 



(4) 



(5) 



is the Hamiltonian of the isolated conductor and finally 

N 



(6) 



X a—1 m—1 
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describes the coupling of the leads to the conductor. The 
set {jam)} represents a basis of scattering states in the 
isolated leads at the Fermi energy Ep- The Hilbert space 
of the cavity is spanned by M localized states \x). These 
two sets of states form a complete basis of the Hilbert 
space of the total system. 

The on-shell scattering-matrix for this system at the 
energy Ep is given by 

S{Ef) ^l~27TiW''GW , (7) 

with the Greens function 

G ^ {Ep - H + iTrWW'')-K (8) 

The matrix elements of S can be written as 

Sam,/3n = SapSmn - 27ri W"J^„GW/3„ , (9) 

where we introduced partial coupling matrices 

Warn = ^ \x){am\Wx,ani ■ (10) 

X 

These matrices describe the coupling of a single channel 
\am) to the cavity. With this definition we can decom- 
pose the M y. M matrix WW^ into a sum 



N 



N 7V„ 



a — 1 a — 1 m— 1 

For later use, we define 

r„ = WaWl and 



(11) 



We want to express the LPDOS in terms of the Hamilto- 
nian Eq. (^). This is possible with the help of expressions 
which relate the LPDOS to functional derivatives of the 
scattering matrix p0|-p2[. The LPDOS of carriers which 
are injected through contact /?, reach point x and are 
emitted into contact a is given by 



u{a, X, 13) 



1 

47ri 



-Tr 



'a/3 



5U{x) 



h.c. 



(12) 



where s^fs is that submatrix of the full scattering matrix 
which describes scattering from all the channels in con- 
tact (3 into all the channels in contact a. The injectivity 
of contact /3 into point x is the sum of all LPDOS over 
all contacts through which a carrier can possibly exit the 
sample, 



1 



z/(x, P) = — ; — : y Tr 



6U{x) 



(13) 



The emissivity into contact a of point x is the sum of 
all LPDOS over all contacts through which a carrier can 
possibly enter the sample, 
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v(a,x) = —- — > Tr 



'a/3 



fea/3 

SU{x) 



(14) 



Finally, the local density of states is given by the sum 
over all injectivities, or a sum over all emissivities, or the 
sum of all LPDOS, 

To find the functional derivation of S with respect to 
U{x), we notice that in the discretized Hamiltonian 
the potential U{x) appears only in the diagonal terms 
Hxx = Ex + U (x) . Therefore we can express the func- 
tional derivative with respect to the potential as an ordi- 
nary derivative with respect to the diagonal elements of 
the Hamiltonian 1261, 



Ssa 



/3 



Ss 



a/3 



SU{x) 



5Hx 



(16) 



Using Eq. (H) for the S'-matrix elements and Eq. (|T^ we 
find 



v{a,x,l3)^Re [5^13 [GFaG]^^ 
-I- 2TTi [GF/jG^FaG]^ 



(17) 



Here, \A\xx denotes the diagonal element of the M x M 
matrix A. Taking into account that the S'-matrix is uni- 
tary we get for the injectivity 



Ka;,/3)=EKa,x,/3)= [GF^Gt]^^ 



(18) 



= / G{x,xi)Tf}{xi,X2)G^ {x2,x)dxidx2 (19) 



(20) 



and for the emissivity 

i^{a, a;) = E iy{a, x, (3) = [G^F^G] 



Note that the injectivity depends in an explicit manner 
only on the coupling elements Wp to the leads through 
which carriers enter and the emissivity depends only on 
the coupling elements to the lead to which carriers finally 
exit. 

The M X M matrices 



= GF^Gt , 
lV„ = G^F^G 



(21) 
(22) 



whose diagonal elements are the injectivity (emissivity) 
are called injectivity (emissivity) operator. However, the 
non-diagonal terms of these operators are not LDOS but 
play an important role in the description of extended 
tunneling contacts. The total LDOS is given by 



j,{x) = J2 ^("' f^) ^ [^^rGt]^ 



(23) 



13a 
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For later reference we also need the probability for trans- 
mission of a carrier from contact a into a different contact 
(3. It is obtained by summing the squared value of the 5*- 
matrix element |s^„_ctmp over all channels m of contact 
a and all channels n of contact /3, 



(24) 



Eq. ( pil) is a standard result that is discussed in text- 
books |M • 



III. INJECTIVITY AND EMISSIVITY AND 
SCATTERING STATES 

In this section we discuss briefly expressions for the 
injectivities and emissivities in terms of the scattering 
states ipam (x) for a system described in the previous sec- 
tion. A scattering state is obtained from a spatially very 
wide (energetically very narrow) wave packet which is 
incident in contact a in channel m in the limit that the 
energy spread tends to zero. The scattering state con- 
sists thus of an incident wave in channel m of lead a and 
typically of waves that are reflected back into all chan- 
nels of lead a and of waves that are transmitted into 
all channels of all the other leads We assume that 
the incident wave is normalized to unity. The scattering 
states of all channels and all leads together with possi- 
ble bound (localized) states form a complete orthonormal 
set of states. To establish a relation between the injec- 
tivity and the scattering states a connection between the 
functional derivatives of the scattering matrix and the 
scattering states is needed. Such a connection was found 
in the analysis of tunneling times based on the local Lar- 
mor clock ||28| . For a more recent discussion we refer the 
reader to Refs. ||2^,^. Here we state only the results. 

The injectivity of a contact (3 into a point x is that 
part of the LDOS which is contributed by the scattering 
states describing particles incident from contact f3, 



1 



^ hvpm{B) 



(25) 



where = ■\/2/m(i? — Eo,pm) is the velocity of elec- 
trons with energy E in channel m of lead (3 and -Eo,/3m 
is the band offset. B is the magnetic field. Using the 
symmetry relation (H) the emissivity can be written in 
terms of scattering states in the form 



v{B; a, x) 



E 



1 



hVa,n{-B) 



\lparni-B;x)\" 



(26) 



With the help of Eqs. ( |25D and (26) we can describe the 
injectivity and emissivity on atomic length scales. The 
spatial variation contained in these densities of states is 



extremely complex not only inside a conductor but even 
in a perfect conductor leading up to a barrier. Consider 
a scatterer characterized by reflection amplitudes ru^rnn 
for carriers incident in contact 1 in channel n and re- 
flected into contact 1 into channel m. In the lead connect- 
ing the scatterer to contact 1 the absolute squared value 
of the scattering state -01, n has then "diagonal" contri- 
butions proportional to Rii,mnXmiy) which are spatially 
independent. Here Rn,rnn = kii, mnP is the reflection 
probability and Xm iv) is the transverse wave function in 
channel m in lead 1. In addition to these diagonal terms 
which are independent of x there are interference terms 
proportional to r*j^^inrii,7nnXi{y)Xm{y)exp{i{ki - km)x) 
which oscillate on a very large length scale for subbands 
with Fermi wave vectors which differ very little. Suppose 
that we are not interested in the very detailed structure of 
the injectivity but only in the injectivity averaged over 
the cross section of the conductor. Integration over y 
eliminates the long range oscillations due to the orthogo- 
nality of the transverse wave functions. The only remain- 
ing oscillations along x are then Friedel-like, proportional 
to cos[2kmx)^ and the longest period is half a Fermi wave- 
length of the topmost occupied subband. Suppose that in 
addition to integrating over the transverse cross section 
we also average the density over a length large compared 
to this period. The resulting injectivity of contact 1 to 
the left of the scatterer is then 



i=l 



i=l 



hvil 



(27) 



and to the right of the scatterer is 



(21) 



•i2 



(28) 



Here the brackets <> indicate the integration over x and 
y. The probabilities i?-"^ = ^n.ij = J2j kn.iiP 



(21) 



Z]j72i,ij = J2j l*2i,ijP are the total probabili- 
ties of all full incident channels which contribute to reflec- 
tion in channel i and to transmission in channel i. Eqs. 
( p7| ) and ( p8[ ) will be useful to discuss the connection of 
the present work with the results of Azbel [|| and Refs. 
pO| , and ||T^ . We only note already at this point that the 
injectivities and emissivities averaged in this way cannot 
exhibit a full Hall effect since the Hall effect depends sen- 
sitively on the variation of the injectivity and emissivity 
in the y direction. 

We have now given expressions for the injectivities and 
emissivities in terms of functional derivatives of the scat- 
tering matrices, in terms of the diagonal elements of the 
injectivity and emissivity operators, and in terms of scat- 
tering states. Our next task is to express the transmission 
probabilities in terms of these densities of states. 
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IV. WEAK COUPLING OF TWO CONDUCTORS 

Consider now two conductors as shown in Fig. ^j. Each 
of the conductors is described by a Hamiltonian Ha and 
Hb, defined according to Eq. (||). The two conductors 
are coupled and this coupling is described by a matrix 



AB 



(29) 



We are interested in the weak coupling limit and, there- 
fore, take the matrix elements I^a^b to be small. The 
Hamiltonian for the total system reads 



"H = TiA + Hi 



AB 



^Ib). 



(30) 



Let Ga and Gb be the Greens functions defined in Eq. 

of the uncoupled systems A and B. Then, to the low- 
est order in ||7^_b|| we can write for the Greens function 
of the coupled system, 

G = Ga + Gb — GaTabGb - GbT^^Ga + 0{\\Tab\\^) ■ 

(31) 

We put this expression into Eq. ( 24 ) to find the transmis- 
sion probability of a contact a of system A to a contact 
/3 of system B, 



T, 



13a 



TabG\t pG bT^qG A^ aG\ 



AB\ 



(32) 



At first glance this formula looks quite cornplicated, but 
if one compares it with Eqs. (21) and (p2[), one sees 
that Tafs is just a combination of the injectivity- and 
emissivity-operators and the coupling matrix Tab , 



0{\\Tab\\') 



(33) 



Let us now assume that the coupling of the two conduc- 
tors is point like. Then we can describe the coupling 
with a single weak bond. Thus, we set Tab = In- 
putting this coupling matrix into Eq. (32), the transmis- 
sion probability reduces to the simple expression given by 
Eq. (0). The probabihty of a carrier to go from contact 
a via the weak link to contact (3 is given by the product 
of the injectivity of contact a into the connecting point 
Xq multiplied by the emissivity of point Xq into contact 
/3. 

Suppose now that one of the conductors, for exam- 
ple conductor B, has only one contact. Then conductor 
B provides a simple description of the tip of an STM. 
(Even though the current distribution between an STM 
tip and the surface is spatially somewhat extended [ p9[ , 
the coupling between tip and surface is theoretically most 



often treated as being point-like.) If there is only a sin- 
gle contact the injectivity and the local density of states 
are identical, UBi^^ tP) = vb{xq). Thus the probability 
for transmission from a contact a of the sample into the 
tunneling tip is given by 



T^a = 4^2^s(xf)|i|2j.^(4,a). 



(34) 



While on the tip-side only the local density of states en- 
ters, on the sample side the relevant density of states is 
the injectivity. If the tip acts not as a carrier sink but 
as a carrier source the transmission probability from the 
tip into the sample is given by 



T^p = An''yA{a,x^)\t\^yB{xE) 



(35) 



and contains on the sample side the emissivity as the 
relevant density of states. 

We see that the transmission is not proportional to the 
total LDOS at the coupling point in the sample, but only 
to the injectivity (emissivity) of that contact a, for which 
we want to know the transmission probability into the 
tip. This is due to the fact that the sample is connected to 
more than one reservoir. If the sample is only connected 
to one electron reservoir, all L(P)DOS of the sample are 
identical and Eq. (|^) gives the Bardeen formula. The 
transmission probability is T = ATT^VA[xQ )\t\^UA{xQ). 
Thus, Eq. (|^) can be seen as a generalization of the 
Bardeen formula for transmission between two multiter- 
minal conductors. Furthermore, Eq. ( p2| ) is the general- 
ization of Eq. (|^) for the case when the tunneling contact 
is not pointlike, but allows for multiple tunneling paths. 
An equivalent expression for transmission between two 
one-terminal conductors coupled via an extended tun- 
neling contact has been given by Pendry et al. |3C|]. 



V. THE VOLTAGE MEASUREMENT 

Now we come back to the case, where the scanning 
tunneling microscope is used to scan along a mesoscopic 
wire and to measure the voltage at different points along 
the wire, c. f. Fig. |l|. The electrochemical potential which 
is applied to the tip reservoir is such that there is no net 
current flowing through the tip. At zero temperature, in 
terms of transmission probabilities, the measured elec- 
trochemical potential is given by |l9| 



M3 



T31 + T 



32 



(36) 



in linear response to the applied potentials ^1 and /i2. 
The transmission probabilities are evaluated at the Fermi 
energy. Tap = Tap{Ep). Putting our expressions for 
the transmission probabilities, Eq. ([l]), into this formula 



gives 



M3 



iy{x) 



^J■2 



^{x, 1) 
v{x) 



(37) 
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First, we remark that since the L(P)DOS depend on the 
position X in the wire where the tip is placed, also the 
measured potential ^3 depends on this position. Second, 
the measured potential does neither depend on the DOS 
in the tip, nor on the coupling strength t between tip 
and sample. All terms in the nominator as well as in the 
denominator arc proportional to the coupling constant t 
and the density in the tip i^Bix) so that these terms drop 
out. (For non-zero temperature, if the DOS in the tip de- 
pends significantly on energy |3f[| , the measured voltage 
depends also on the density of states of the tip). 

This formula allows us now to assign to every point on 
the wire an electrochemical potential. However, there is 
no simple relation between this measured electrochemi- 
cal potential and the electrostatic potential in the wire. 
The electrostatic potential itself can not be measured, 
at least not using the method described here. Measuring 
the electrochemical potential /i3 at a certain point x does 
not mean that the electrons at that point x locally are 
distributed according to a Fermi function with the elec- 
trochemical potential fi^. We emphasize that there is no 
inelastic scattering inside the sample. Applying a bias 
M2 ~ Ml brings the system to a non-equilibrium state so 
that the electrons inside the sample are not distributed 
according to a Fermi distribution. Next we illustrate the 
content of this formula with two examples. 



A. Friedel-like oscillations across an impurity 

Eq. ( ^tI) is valid for any distribution of impurities in 
the wire. The only problem is to find expressions for the 
L(P)DOS. For a complicated geometry with randomly 
distributed impurities there is no hope to find an exact 
analytic expression for the L(P)DOS. But we can con- 
sider a simplified example, which can give an idea how 
the measured potential should look like in the neighbor- 
hood of an impurity. 

We consider a one channel perfect conductor which 
contains only one scatterer. Discussions closely related 
to the point of view taken here are given by Levinson | |3^ ] 
and in Ref . . We assume that the equilibrium poten- 
tial is constant all along the wire except for a delta-peak 
at X = 0. On the left side, the wire is connected to elec- 
tron reservoir 1 , and to the right side to electron reservoir 
2. For this model |2|], it is easy to find the analytic form 
of i^ix) and i'{x, 1). In Fig. ^ both quantities as well 
as their ratio are shown as functions of the position x. 
The LDOS, which consists of contributions of scattering 
states coming in from the left and from the right side, 
oscillates on both sides of the scatterer. The injectivity 
of contact one consisting only of contributions of scatter- 
ing states coming from the left side oscillates only to the 
left of the scatterer. Only there we have interference of 
incoming and reflected waves. To the right of the scat- 
terer, we have only outgoing waves so that the injectivity 
of contact one is constant. 



2 - 



1 - 




-1 1 

m 

FIG. 3. The LDOS u{x) (solid line), injectivity iy{x,l) of 
contact 1 (dotted line) (both in units of 1/hv), and the ra- 
tio of these two densities vix, l)/i'{x) (dashed line) of a one 
dimensional wire with a S barrier at a; = which leads to a 
transmission probability of T = 0.8. 

The X dependence of the measured potential comes 
from the oscillations in the ratio of the two densities, Eq. 
(^). Since the injectivity is a part of the LDOS, this 
ratio is always between zero and one and thus, the mea- 
sured potential 113, Eq. lies always between /zi and 
fi2- If the ratio is close to one, which is often the case 
to the left of the scatterer, the measured potential /13 is 
close to the applied electrochemical potential to the left 
of the scatterer ^1. However, we can also find positions 
to the left of the scatterer where the ratio is small so that 
the measured potential is close to the electrochemical po- 
tential to the right of the scatterer ^2- Likewise, we can 
also find positions to the right of the scatterer where we 
measure a potential which is close to the applied poten- 
tial at the left side. This leads to interesting effects when 
the voltage probe is used to measure the resistance of a 
barrier. 



VI. THE RESISTANCE MEASUREMENT 



scatterer 



|Lti 





T 


X4 r 











|Lt2 



1^3 ' ILL4 

FIG. 4. Experimental setup for the measurement of the 
four-terminal resistance of a scatterer using an STM tip as a 
voltage probe. 

The four-terminal resistance (Fig. ^ of a scatterer 
is defined as the ratio of the voltage drop across the scat- 
terer, /i3 — /J.4, divided by the current flowing through the 
scatterer due to the applied voltage /ii — /i2, 



7^ 



12,34 



M3 - M4 
e/12 



(38) 
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h 1 (731X42 - 732^41) 

e2r(T3i+T32)(T4i+T42) 



(39) 



Here T — T7-(si2s|2) is the overall transmission proba- 
bility of the conductor in the absence of the measuring 
contacts. Using Eq. (|l|) the four terminal resistance is 



12,34 



h 1 /(1/(X3,1)1/(X4,2) - Z^(X3,2)Z^(X4,1) 



T 



h 1 f j/(x3, 1) 



;/(x3)z^(x4) 

1'(X4, 1)' 



T \ I/(X3) J/(X4) 



(40) 
(41) 



Notice that the last expression is just the difference of two 
(three terminal) voltage measurements given by Eq. ( 37 ) . 
The notion x^ = {xi,yi),i = 3,4 denotes the coupling 
points of the voltage probes 3 and 4 in the x — y plane of 
the wire. Since in the densities j^(x, a) and i^(x) interfer- 
ence between incoming and reflected waves is taken into 
account, Eq. ( ^2| ) is a phase-sensitive resistance. Due to 
these interference effects the densities show a complicated 
spatial behavior, e. g. the Friedel-like oscillations in the 
one-channel case. 

Instead of using the weak coupling contacts as voltage 
contacts we might also use contact 3 to inject current and 
contact four as the current sink. The measured resistance 
is then 7^34 12 and is related to 7^12,34 by a reciprocity 
relation 7^34,l2(-B) = 7ei2,34(-'B). Thus this resis- 

tance is determined by the difference of the emissivities 
into contact 1 of the points X3 and X4, 



7?-34,12 - 



h 1 f t^(l,X3) Z^(1,X4) 



J/(X3) 1/(X4) 



(42) 



Note that both the overall transmission probability T 
and the local densities i^(x) are even functions of the 
magnetic field. We also remark that one might believe 
that current transport from one weak coupling probe to 
another invokes more information then is contained in 
the injectivities or emissivities. This not the case, since 
the current balance is such that to lowest order in the 
coupling strength |ip of the weak coupling contacts the 
injected current first reaches the massive contacts 1 and 
2 and the current at probe 4 is determined by carriers 
injected by the massive contacts 1 and 2. Direct trans- 
mission of carriers from one weakly coupled contact to 
the other one is a second order effect, proportional to 
and thus, generally only a small perturbation. Chan and 
Heller showed recently ^] that, on surfaces with point 
defects (adatoms), even this second order effect can be 
deduced from single tip measurements. 

In order to get from Eq. (|4|) to a position indepen- 
dent value for the four-terminal resistance we average 
the phase-sensitive result by moving the voltage probes 
on both sides of the scatterer over some distance while 
measuring the voltage. One possibility is to keep the 
transverse coordinates j/q fixed and average only along 
the X axis. 

Let us now consider the one-channel case. There we 
know that it is sufficient to average over half a Fermi 



wavelength, since the measured voltages show a peri- 
odic behaviour on this length scale. We get the phase- 
averaged result. 



.hi 
('7^12,34) — —77; 



vjxs., 1) 
1^(2^3) 



yjxj, 1) 
v{xa) 



h 1-Vt 
? f 



(43) 
(44) 



The same result was already found by Biittiker who 
described the voltage probes as wave splitters. 

For comparison we calculate also a phase-insensitive 
resistance which means that we neglect the phase coher- 
ence of incoming and reflected wave altogether. This is 
equivalent to averaging injectivity and local density sep- 
arately and leads to the Landauer formula for the 
resistance of a scatterer, 



TZ 



h 1 f{iy{x3,l)) {iy{x4,l)) 



12,34 



(K^3)) 



{iy{x4)) 



hl-T 



(45) 
(46) 



We emphasize that what can be measured directly by us- 
ing a sufficiently sharp tip is the phase-sensitive result. 
By moving the tip and averaging over the measured po- 
tentials we get the phase-averaged result. Note that the 
phase-averaged result would also be obtained if the mea- 
surement is made further then a phase-breaking length 
away from the scatterer. 



A. The few-channel case 



For a scatterer connected to leads with iV > 1 open 
channels Eq. (^) is still valid. The injectivity as well 
as the local density are now not anymore periodic func- 
tions. They consist of a superposition of oscillations with 
different wavelengths. In fact, the functions contain oscil- 
lations with wave vectors k = ki± kj given by all possible 
combinations of the Fermi wave vectors ki of the N chan- 
nels. If the number of channels is very large, one expects 
that the densities become nearly constant as a function 
of the position x. Instead of the exact densities we can 
then use the densities averaged over a portion of the con- 
ductor. The injectivity of contact 1 to the left and to 
the right of the scatterer averaged over x and y are given 
by Eqs. ( p7| ) and (p8|). Using these densities in Eq. ( |45| ) 
gives the result of Azbel |^ and Biittiker et al. found 
with the help of a charge neutrality argument. 

Since the densities averaged over the x and y coordi- 
nates do not anymore exhibit a dependence on the trans- 
verse y coordinate, the resulting resistance formula can 
not explain the magnetic field dependence (Hall effect) of 
the measured resistance. In contrast, Eq. (p2|), includes 
the exact, spatial densities and shows a dependence on 
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the magnetic field not only through transmission prob- 
abilities but also due to the magnetic field dependence 
of the injectivity and emissivity. We will now illustrate 
this point by investigating the Hall resistance of a one 
channel wire with an obstacle. 



B. Magnetic field dependence of tlie Resistance 



Let us consider a scatterer which is connected via ideal 
leads to two electron reservoirs. Let us assume that 
in the ideal lead, far away from the scatterer, we have 
a uniform potential in the longitudinal x direction and 
a parabolic confining potential in the transverse y di- 
rection, U{y) = l/2mcjQy^. Furthermore, the lead is 
threaded by a magnetic field perpendicular to the x — y 
plane. In the lead the eigenfunctions of such a system 
can be written as a product of a plane wave cyi'p{ikx) in 
X direction and a transverse wave function Xk{y)- The 
transverse wavefunction depends now on the k vector of 
the plane wave in the x direction. That means, that not 
only different channels, but also incoming and refiected 
waves in the same channel have different transverse wave 
functions. This mechanism leads to a spatial separation 
of incoming and reflected waves and in a strong magnetic 
field to the formation of edge channels |^,^. The Hall 
effect in perfect ballistic wires has been discussed H^] 
in connection with the suppression of the Hall effect in 
ballistic crosses jSTj . This suppression is, however, an ef- 
fect which depends on the geometry of the cross p8| , p9[ . 
Here, the main effect which we investigate arises due to 
the scatterer in the wire which is taken to have a mag- 
netic field independent transmission probability T — 0.5. 

For the case of only one open channel in the lead, the 
scattering state coming in from contact 1 can be written 
in the lead connecting the scatterer to contact 1 as 



^i(x) 



(y)+rne-^*-"X-(y) 



(47) 



Here, rn is the reflection amplitude for reflection from 
contact 1 into itself, fc+ and fc_ are the wave vectors for 
incoming and reflected waves and x+ ^-nd X- ^^^^ the cor- 
responding transverse wave functions p7[ |. The voltage 
fi3 measured on a point to the left of the scatterer is given 
by Eq. (|^. Using Eq. we find for the densities, 

Kx,l) = -^(|x+(y)P + |rnnX-(y)P 
nv 

+2\rnX+{y)X-{y)\cos{2kx + S)) , (48) 

Kx,2) = -l(l-|rnP)|X-(y)P, (49) 
nv 

i^(x) = z/(x, 1) -f i/(x, 2) . (50) 

Here, k = (fc_|_ + k-)/2 and rn = |rii| exp(zi5). The mag- 
netic field dependence of this formula is hidden in the 
transverse wavefunctions. 
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FIG. 5. Magnetic field dependence of the resistance of a 
ballistic wire with a scatterer with transmission probability 
T — 0.5. The voltage probes are placed on opposite edges 
of the lead having the same x coordinate (solid line), being 
separated in x direction by a quarter of a Fermi wavelength 
(dotted line) and (dashed line). 

In Fig. I we show the measured four-terminal resis- 
tance as a function of the ratio tOc/i-^o, where lUc ~ eB /m 
is the cyclotron frequency and B the magnetic field. The 
three curves correspond to three different configurations 
of the voltage probes, contact 3 and 4. The two voltage 
probes are placed on the two edges of the lead connecting 
contact one to the scatterer. First in a way that the x po- 
sition of the two probes is the same (solid line). Second, 
the X coordinate of probe 3 is such that the LDOS has a 
minimum, cos(2fca; + 5) = —1, at that point in the lead, 
whereas probe 4 is place on a point where the LDOS has 
a maximum, cos(2A:a:-|-i5) = 1, (dashed line). And finally, 
probe 4 is placed over a minimum of the LDOS and probe 
3 is placed over a maximum of the LDOS (dotted line). 
Note, that, if there is no magnetic field present, the resis- 
tance is zero, when the positions of the two probes differ 
only in their y coordinate. If the magnetic field is turned 
on, a Hall-resistance develops which at strong magnetic 
fields reaches the quantized value TZ ~ h/e^. In the case, 
where the positions of the probes differ in the x coordi- 
nate in the way described above, a non-vanishing resis- 
tance is already measured at zero magnetic field. How- 
ever, when the field is turned on, both curves approach 
the value h/e^. For small magnetic fields, we can expand 
the resistance formula, Eq. ( |4l| ) , and get to first order in 



7^l2 



h 1 ( V(3{xi,2)\ vo(x^,2) 



y4 t^o(a:4,2) ^ 2/3 i^o{x3,2) 
^^0 \yo i^o{xa) yo fo{x3) 



(51) 
(52) 



where lOc — eB/m is the cyclotron frequency and yo — 
y^h/ (mwo)- In this formula, the y dependent Hall resis- 
tance and the x dependent, longitudinal resistance enter 
additively. 
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VII. DISCUSSION 

In this work we described systems consisting of two 
weakly-coupled multi-probe conductors starting from the 
global scattering matrix of the whole system which covers 
all parts including the weakly coupled contact. We de- 
rived a general transmission formula, Eq. (|3^), for trans- 
mission through the weak-coupling contact. We have in- 
vestigated this formula in the case where there is only one 
tunneling path. Using these expressions we can rewrite 
formulae, which express the resistance as functions of 
transmission probabilities, as functions of the local par- 
tial density of states. Applying this result the resistance 
of a one-channel conductor with a barrier shows inter- 
esting and surprising features. Our resistance formula 
can also account for the Hall-resistance. This point is 
illustrated using two weakly-coupled voltage-probes to 
measure the magnetic-field dependence of the resistance 
of a ballistic one-channel conductor with a barrier. 

Based on the general expression, we can also treat 
contacts which permit multiple tunneling paths. After 
all, even an STM exhibits a current distribution with a 
certain spatial width |^^. It is then interesting to ask 
what the densities are which are measured by spatially 
extended tunneling contacts. To our knowledge, a de- 
tailed study of extended tunneling contacts has not yet 
been done. 

We have treated the zero-temperature limit. At a non- 
vanishing temperature the corresponding results are ob- 
tained by multiplying the transmission probability with 
the Fermi function of the injecting reservoir. Even in 
the limit of a single tunneling path the corresponding re- 
sistances will then in general not be independent of the 
density of states in the tip if this density exhibits a sub- 
stantial variation at the Fermi energy ||3l|| . 

We have already emphasized that the density of states 
discussed here are quantities which are conjugate to the 
electrochemical potential of a contact . At zero tem- 
perature they are evaluated at the Fermi energy in the 
equilibrium potential. The densities used here are thus 
essentially chemical quantities and interactions enter only 
through the equilibrium potential. The interaction in- 
duced portion of the density does not enter into the 
transmission behaviour of a conductor. This should be 
contrasted with the notion of tunneling-density of states 
which are evaluated using the unrestricted Green's func- 
tions containing the full interaction. An investigation of 
this important point is beyond the scope of this work. 
We refer an interested reader to a discussion of the same 
issue concerning not density of states but directly the 
conductance of interacting systems fiot] . 

We have discussed examples of voltage and resistance 
measurements on conductors with a single open channel. 
The approach discussed here, however, is also suitable for 
conductors having several or even a large number of open 
channels. For a metallic diffusive conductor of length L, 
extending from x = to x = L, with a local density of 



states v the injectivity can be separated into an ensemble 
averaged part given by i/(l — x/L) and a fluctuating part. 
The average behavior gives us the linear voltage drop and 
the ohmic length dependence which we expect for such 
conductors. More interesting is the investigation of the 
Hall conductance: although the local, ensemble averaged 
density of states of a metallic conductor is independent 
of the magnetic field, the injectivities exhibit a linear de- 
pendence on the magnetic field and and this gives rise to 
the Hall resistance, similar to the one-dimensional exam- 
ple discussed in Section VI. 

To conclude we emphasize that the investigation of 
tunneling contacts on mesoscopic conductors is an in- 
teresting subject which has so far found only limited at- 
tention. 

This work was supported by the Swiss National Science 
Foundation. 
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